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We analyze two approaches to quantum state transfer in solid-state spin systems. First, we 
consider unpolarized spin-chains and extend previous analysis to various experimentally relevant 
imperfections, including quenched disorder, dynamical decoherence, and uncompensated long range 
coupling. In finite-length chains, the interplay between disorder-induced localization and decoher- 
ence yields a natural optimal channel fidelity, which we calculate. Long-range dipolar couplings 
induce a finite intrinsic lifetime for the mediating eigenmode; extensive numerical simulations of 
dipolar chains of lengths up to L = 12 show remarkably high fidelity despite these decay processes. 
We further consider the extension of the protocol to bosonic systems of coupled oscillators. Second, 
we introduce a quantum mirror based architecture for universal quantum computing which exploits 
all of the spins in the system as potential qubits. While this dramatically increases the number 
of qubits available, the composite operations required to manipulate "dark" spin qubits signifi- 
cantly raise the error threshold for robust operation. Finally, as an example, we demonstrate that 
eigenmode-mediated state transfer can enable robust long-range logic between spatially separated 
Nitrogen- Vacancy registers in diamond; numerical simulations confirm that high fidelity gates are 
achievable even in the presence of moderate disorder. 

PACS numbers: 03.67.Lx, 03.67.Hk, 05.50.+q, 75.10.Dg 



INTRODUCTION 



The ability to perform quantum logic between remote 
registers has emerged as a key challenge in the quest for 
scalable quantum architectures HUH- Qubits, the fun- 
damental building blocks of such an architecture are of- 
ten benchmarked by their coherence times [5H7]. Nat- 
urally, those qubit implementations which possess the 
longest coherence times also interact most weakly with 
their local environment, making multi-qubit quantum 
logic in such systems difficult [3j [8]. As a result, there 
has been tremendous recent interest in quantum data 
buses, which enable universal gates between physically 
separated quantum registers [9HT6]. Such data buses 
have been proposed in systems ranging from trapped ions 
[17-19 and superconducting flux qubits [20-22 to cou- 
pled cavity arrays [23H25] and solid-state spin chains [26]- 
[37] . Prior proposals have focused on achieving perfect 
state transfer using either initialized [26, 28, 29, 37 , en- 
gineered [31, 38, 39 or dynamically controlled quantum 
channels 



By contrast, here, we analyze a general method 
for high-fidelity quantum state transfer (QST) using 
an infinite-temperature (unpolarized) data bus [3] [27] . 
Our method requires neither external modulation dur- 
ing state transfer, nor precisely engineered coupling 
strengths within the bus, making it an ideal candidate for 
solid-state spin-based quantum computing architectures 
[U [45] . We envision the long-range coherent interaction 
between remote qubits to be mediated by a specific col- 



lective eigenmode of the intermediate quantum data bus. 
In the solid-state, such eigenmodes naturally suffer from 
localization effects associated with lattice imperfections 
and disorder [271 EH]- Exploration of the interplay be- 
tween such localization effects and intrinsic constraints 
set by finite coherence times, is important to assess the 
feasibility of proposed architectures. 

Our paper is organized as follows. In Sec.[TTJ we extend 
the previously proposed notion of eigenmode-mediated 
quantum state transfer [27] to the transverse field Ising 
model. In addition to being closely related to the ac- 
tual achievable Hamiltonian of certain driven spin sys- 
tems, this simple model enables an analytic description 
of the state transfer protocol. In Sec. |III[ we build upon 
these protocols and derive analytic expressions charac- 
terizing the channel fidelity for state transfer between re- 
mote quantum registers. Next, we generalize our method 
to bosonic systems (e.g . coupled cavities and pendu- 
lum arrays) in Sec. IV In particular, we demonstrate 



that chains of coupled oscillators can faithfully trans- 
port quantum information even at high oscillator tem- 
perature. 

Having described eigenmode-mediated QST in both 
the fermionic and bosonic context, we then turn to a spe- 
cific implementation within a solid-state quantum com- 
puting architecture. In Sec. [V] we analyze eigenmode- 
mediated quantum state transfer between remote spin- 
quantum-registers [47-49 . To be specific, we consider 
Nitrogen- Vacancy (NV) defect center registers and ex- 
amine the optimization of state transfer fidelities in the 
presence of both disorder and a finite depolarization 
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(Ti) time. The interplay between disorder- induced local- 
ization and decoherence yields a natural optimal chan- 
nel fidelity, which we calculate. Ultimately, this opti- 
mization demonstrates the feasibility of scalable archi- 
tectures whose remote logic gates can harbor infidelities 
below the threshold for error correcting codes [5(3 I5T] . 
While eigenmode-mediated QST fundamentally requires 
the register-chain coupling to be weaker than the intra- 
chain coupling, we demonstrate in Sec. |VI| that general- 
izations to the strong coupling regime are also possible. 
We provide numerical simulations in parallel with the 
analytic channel fidelities derived in Sec. |III| 

In Sec. |VII[ we perform exact diagonalization for spin 
systems, which includes the full long-range dipolar inter- 
action. We find remarkably high fidelities for our pro- 
posed QST protocols in chains of length up to L = 12. 



Finally, in Sec. VIII, we describe and analyze an alternate 
architecture, which utilizes globally controlled pulses for 
state transfer [40j[41]. In this case, we demonstrate that 
all spins in the system (e.g. even dark intermediate chain 
spins) can be viewed as potential qubits. However, while 
this dramatically increases the number of qubits avail- 
able, the composite operations required to manipulate 
such intermediary spin qubits significantly raise the er- 
ror threshold for robust operation. 



II. EIGENMODE-MEDIATED QST 

In this section, we begin with an idealized system 
in which to understand eigenmode-mediated QST [27] , 
namely, the transverse field Ising model, 



H 



N-l 



N 



(i) 



where n is the nearest-neighbor coupling strength and 
B represents a uniform transverse field on each site. In 
addition to being realizable in a variety of experimen- 
tal systems, ranging from NVs and trapped ions to elec- 
trons floating on helium [27j [52j [53] , this model also has 
the virtue of being exactly solvable; this will allow us to 
clearly illustrate the essence of eigenmode-mediated state 
transfer and to understand the many-body entanglement 
which arises. 

Expanding <jf as a function of spin flip operators, 
erf = (erf ± zcr^)/2, and utilizing the Jordan- Wigner 



fermionized Hamiltonian, 



N-l 



, yields the 
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which is quadratic and conserves fermionic parity with- 
out conserving particle number. To solve Hjw, 



we re-express it as 
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via an orthogonal matrix, O, such that OAO T = A. The 
eigenmodes come in pairs with energy dze^, corresponding 
to eigenvectors dk = 2 k-ij4>j an d d\ = 2 k,j4 > ji where 
k = 1, • • • , N. This transformation yields 



N 



H JW = ^2 e k( d k d k ~ d kd\), 



(4) 



k=l 



where the d-modes satisfy standard Dirac anticommu- 
tation relations. For a uniform chain the spectrum is, 
€k ~ y/ft 2 + B 2 — 2Bk cos where — kir/(N + 1). 

We now consider the addition of quantum registers, 
labeled and N + 1, at the ends of the data bus (Fig. 
[I]). The registers couple perturbatively with strength g 
to the ends of the Ising spin chain [27] [28] and we apply 
a local Zeeman field B' , 

H' = -g(o$of + a%a% +1 ) + B'(a z + a* N+1 ). (5) 



Upon fermionizing, 
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By tuning B' = e z , we ensure that the external regis- 
ters are coupled resonantly to a single finite-energy eigen- 
mode d\ of the intermediate chain. Quantum state trans- 
fer proceeds via resonant tunneling through this mode. 

Noting that a = Ef=i(° T k2/c-i4 + E^=i(° T )i,2/c4 
allows us to re-express c± and cat in terms of the d- modes. 
By choosing g0 2z -i i i = g0 2z -i,N < B' , \e z - e z±1 \ we 
ensure that off-resonant eigenmodes are only weakly cou- 
pled to the quantum registers, leaving an effective three- 
mode picture, 

H e ff ~ e z {d\d z - d z d\) + e z (c\co - c cj) 
+ e 2 (4v+i c Ar+i - cat+i4v+i) 

- g0 2z -i,i(cld z + 4 c o) 

- g0 2 z-i,N(J N+1 d z + dlc N+1 ). (7) 

It is interesting to note that for B < ft, the Hamilto- 
nian in Eq. (2) represents a spin-less p-wave supercon- 
ductor in its topological phase [55]. The zero energy 
boundary modes of this system have received a great 
deal of attention recently. As these Majorana zero modes 
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FIG. 1: (color online). Quantum state transfer is achieved 
by tuning the left and right quantum registers (blue) to a 
particular eigenmode (red) of the intermediate data bus. By 
ensuring that the coupling, t between the registers and the 
chain is sufficiently weak relative to the spacing of adjacent 
eigenmodes, it is possible to consider evolution in an effective 
three-mode picture. Such eigenmode-mediated QST is ap- 
plicable in a variety of contexts, ranging from solid-state spin 
chains to coupled bosonic degrees of freedom (e.g. pendulums 
or cavity arrays). 



are exponentially localized, they cannot be employed for 
state transfer. In our analysis, this follows from the fail- 
ure of the secular approximation to remove fermion num- 
ber non-conserving terms. A straight-forward calculation 
shows that the pairing terms precisely cancel the hopping 
terms in the effective evolution. 

Equation (7) represents the key result of the above ma- 
nipulations. State transfer is achieved by time-evolving 

for r = 7= ^ , leading to unitary evolution, 

v 2go 2z -i,i 

JJ e ff = e~ irHeff = (-l) n -(-l)( c o+ c Wi)( c o+Civ+i)/2 

= (-^-(l-fcJ + c^Ccd + cjv+i)), (8) 

where n z = d\d z . It is instructive to write the ex- 
plicit action of U e ff on the subspace spanned by \I/ = 
cj v+1 |^),CQcj v+1 |^)}, where \Q) is the vac- 
uum associated with cq, c/v+i, 



u eff 9 = (-iy 
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Up to signs, the effective evolution in the register sub- 
space is a swap gate. In the spin representation, ow- 
ing to Wigner strings, there exists an additional set of 
controlled phase (CP) gates, as shown in Fig. [2] Since 
CP 2 = I, this entanglement can be easily cancelled and 
logic gates between the remote registers can be success- 
fully accomplished [27j HSJ [56] . We detail two possible 
such protocols below. 

One protocol, herein termed "there-and-back" , is par- 
ticularly applicable to the case of multi-qubit quantum 
registers. For a two-qubit register, we can label one 
qubit as the memory qubit while the other represents 
the "coupling" qubit. Once an eigenmode-mediated swap 
between the coupling qubits is accomplished, an intra- 
register CP-gate is then performed between the two 



qubits of the remote register. The return swap then 
cancels the unwanted entanglement illustrated in Fig. [2j 
leaving only a controlled-phase gate between the two 
memory qubits. Since CP gates, together with single- 
qubit rotations, can generate arbitrary unitary opera- 
tions, such a procedure enables universal logic between 
remote registers. 

An alternate method, which we call the "paired proto- 
col" utilizes a two-qubit encoding to cancel the Wigner 
strings. In this approach, the quantum information is 
encoded in two spins, a and 6, with logical basis | \) = 
I l)a| 1)6, | t) = I t)a| t>6 [131451156]; the intuition be- 
hind this encoding is that it produces an effective bosonic 
excitation, thereby mitigating the effect of the fermionic 
Wigner strings. State transfer proceeds by successively 
transferring a and b through the intermediary chain. 



III. ANALYTIC CHANNEL FIDELITY 

We now derive the channel fidelity associated with the 
paired protocol. To set up the analytic framework, we 
begin by calculating the fidelity of a simplified protocol, 
termed the "double-swap". In this double- swap, we con- 
sider the left register (indexed 0) undergoing two succes- 
sive eigenmode-mediated swap gates. Ideally, this sim- 
plified protocol swaps the quantum information twice, 
thereby disentangling it from the intermediate chain and 
also returning it to its initial position at the left register. 
We then consider a second protocol, termed the "single- 
swap" , in which the quantum information undergoes only 
one eigenmode-mediated swap-gate. Analyzing this pro- 
tocol will illustrate the effect of the residual entanglement 
on the channel fidelity. Finally, we turn to the paired- 
protocol and demonstrate that the proposed two-qubit 
encoding can eliminate this entanglement, thereby en- 
abling quantum state transfer. In Appendix [B] we com- 
pute the channel fidelity for an eigenmode-mediated re- 
mote a z gate. 



A. Double-swap 

The average channel fidelity for a quantum dynamical 
operation is given by 



1 

12 



E Tr[<7^)] 



(10) 



where £ characterizes the quantum channel [57 . For 
simplicity of notation, we will restrict ourselves to the 
XX-model [27], H = g(cr^a^[ + <7^cr^ +1 + h.c.) + 

J^^ 1 n{crfcr~ +1 +h.c), although analogous results hold 
for the previously considered transverse field Ising model. 
For the double-swap (DS), we let U represent evolution 
under H for a time, t = 2r, equivalent to twice the state- 
transfer time. Let us suppose that the left register is ini- 
tially disentangled from the remainder of the chain, which 
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FIG. 2: (color online). Schematic circuit diagram depict- 
ing eigenmode-mediated state transfer between the quantum 
registers (QR) and N + 1. Controlled phase gates are repre- 
sented as circle-ending dumbbells while X-ending dumbbells 
depict a swap gate. In addition to the desired state trans- 
fer, each register is CP-entangled with all intermediate spins 
owing to the Wigner strings associated with fermionization. 
This additional entanglement can be cancelled by utilizing a 
simple two-qubit encoding. 



is in a thermal mixed state p^f] the average double-swap 
channel fidelity is then given by, 
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E Tr[,j[/K^ c f)[/t] 
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-x,y,z 

= \ + Y2 J2 Tr[4(t)(4®p? h s )], (11) 

i=x,y,z 

where <J l ^{t) is the Heisenberg evolution of the left regis- 
ter. By fermionization, this evolution can be re-expressed 
with respect to elements of the matrix M = e~ lKt where 
K is the (N + 2) x (N + 2) coupling matrix of the full 
Hamiltonian (including registers), H = J2f^ Kijc\cj- 
Evolution of the fermi operators is governed by c m = 
-iJ2n K rnnC n , implying that c m (t) ^ n M mn c n and 
further, that 

cr+(£) = U ] (j+U = U ] clU 

= E M o*4 = E M ^t n e<TO,+ff, ~ > ( 12 ) 

i i l<i 

a z (t) = 2cl(t)c (t) - 1 = -1 + 2 J2 M&Mojclcj 

ij 

= -l + 2^M *M„^: J] ^ + ^,(13) 

ij i<l<3 

where we have used the fact that c\ carries no 
Wigner string. To evaluate Fjjs, we n °te that 
a ± = (a x ±ia y )/2, and hence, Tr [(T%(t)(a% <g> p ch )} = 



Tr [(cr+(t) + a (t))((cr+ + cr Q ) (g) p ch )] . Contributions 
are only obtained from the cross-terms, 0"o~(t)(0o~ <g> p c h) 
and (Jq (t)(crQ pch), since the number of excitations in 
i = must be preserved to generate a non-zero trace. 
For example, using Eq. ( [l2| ), 

Tr [a+(t)(aQ®p ch )] 



Tr 



(E M o>, + ri ( 



)O ®Pc/i) 



= Tr [M >+a - ® Pch ] = M* . (14) 
An analogous calculation yields Tr [a^ - (£) (<7q~ (8) p c h)] 



Moo- Finally, for the a z terms, one finds, using Eq. (13), 
Tr [a z (t)(a z p ch )} = Tr [-a* p ch ] 
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Tr [2M* M 00 a+a -ag p c/l ] = 2|M ( 



oo I 



(15) 



where we've noted that i = j to ensure that the number of 
excitations in each mode is conserved. Moreover, we must 
also have that i = j = 0, since Tr [erg] = 0. Combining 
the above terms yields the double-swap channel fidelity 

as, 
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FdS = 2 + 6 (M °° + M °*° + |M ° o1 )• 



(16) 



Interestingly, we need to compute only a single matrix 
element to obtain the relevant channel fidelity. 



B. Single- swap 

We now consider the single-swap (SS) channel fidelity 
associated with the transfer of quantum information from 
the right register (indexed N + 1) to the left register 
(indexed 0), 

Fss = l + ^ J2 Tv[al(t)(p s J®^ N+1 )}, (17) 



z=x,y,z 



where p^f now characterizes the initial state for spins 
{0, • • • , N}. Note that Fss will be independent of the 
direction of state transfer, and we have chosen right to 
left for notational simplicity. From Eq. (12), one finds, 



°oV) = + c o (*) = E M oi4 + M ^ 
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= E[{MM)*K + lm(M 0i )*y} H(-a?){18) 
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In analogy to the DS case, i ^ N + 1 terms do not con- 
tribute to the trace, 
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HK(t)(Pc* ® a% +1 )] = 2Re(M 0tN+1 )Tr[p s J l[(-af)]. 

(19) 
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The a y term yields an identical contribution while the a z 
term yields, Tr[a$(t)(p^ ®a z N+1 )] = 2|M (V v+i| 2 . There- 
fore, 



Fss = l + l[2Re(M , N+1 )Tv[p s JY[(-af)} + \M , N+1 



N 



2 6 



1=0 



(20) 

For perfect transfer with Fss = 1, we would require 

both \M 0iN+1 \ = 1 and \T*[p S JIlZ (-°t)}\ = 1- In 
the case of an unpolarized chain, the second condition is 
unsatisfied since the expectation value of the chain par- 
ity operator P = n^o( — a i ) ^ s zero - The dependence of 
the single-swap fidelity on the intermediate chain's parity 
demonstrates the entanglement illustrated in Fig. [2j and 
presents an obvious problem for QST. 



Paired-Protocol 



one must take care to correctly evaluate such strings, an 
otherwise straightforward computation yields, 



F enc = -(2|M ,Ar+i| 2 Re [M^ N+1 - M 0i0 M N + liN + 1 ] 

N 



+ |M ,j>r+i| 2 + |^M J >r + i,<M i ,o| 2 ) + 5. (22) 
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Again, one only needs to compute certain matrix ele- 
ments of M, and, in fact, an analytic form for all such 
elements can be obtained (see Ref. [28 and Appendix |A|). 

Before concluding this section, we point out that one 
can alternatively decode the quantum information into 
qubit (TV + l) a via C(7v+i) a NOT(7v+i) 5 - In this case, the 
expression for F enc is identical to Eq. (22), except the 



term | ^f =1 Mn+i^Mi^\ 2 is not present. Thus, the de- 
coding into qubit (TV+1)& described above gives a slightly 
higher average fidelity and we will use this decoding in 
later numerical simulations. 



To overcome this problem, we now turn to the two- 
qubit encoding proposed in Sec. [ll| i.e. | i) = | i 

)a\ i)b, | t) = | t)a| t)&- Let us index the ful1 

chain as {0 a , 6 , 1, • • • , TV, (TV + 1) 6 , (TV + l) a } and de- 
fine Ub as the transfer process through the sub-chain 
{Ofr, 1, • • • , TV, (TV + 1)5}, while U a represents the transfer 
process through the sub-chain {0 a , 1, • • • , TV, (TV + l) a }- 
To model a realistic experimental situation, we will as- 
sume that the quantum information is originally encoded 
in qubit a , while qubit is prepared in state | A 
Co a NOT 0b gate is then applied to encode the informa- 
tion in the logical 0-register. After the state transfer, 
we apply C(Ar+i) b NOT(7v+i) a to decode our quantum in- 
formation into qubit (TV + 1)5. The unitary character- 
izing the encoding, state transfer, and decoding is then 
U=C (Ar+i) b NOT( A r+i) a ^6^aCo a NOT 0b , and the average 
channel fidelity is given by 
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E Tr Niv + i) b w(<^itk(ti 
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Pch ®PN+l) ■ 
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Here, p^ is the mixed initial state of the intermedi- 
ate chain ({l,--- , TV}), while pN+i is the mixed state 
of the encoded (TV + 1) register within the logical sub- 
space. Working within this logical subspace is crucial 
to ensure that CPo a ,jv+i a CPo b ,jv+i b = I. Inspection 
reveals that the paired-protocol includes two contribu- 
tions from the chain parity operator, and since P 2 = I, 
we have effectively disentangled from the intermediate 
chain. Since a consistent ordering of the spin-chain is re- 
quired to implement the Jordan- Wigner transformation, 
the Hamiltonian, Hjj a governing the U a transfer evolu- 
tion will contain uncanceled Wigner strings. For exam- 
ple, the piece of Hu a containing the coupling between the 
registers and the ends of the spin-chain takes the form, 
H Ua =^(4 a e^c 1 +cj v e™+ 1 ).c (iV+1)a +h.c.). While 



IV. GENERALIZATION TO OSCILLATOR 
SYSTEMS 

In this section, we analyze the generalization of 
eigenmode-mediated state transfer to systems of bosonic 
oscillators. The realization of such coupled-oscillators is 
currently being explored in systems such as, cavity arrays 
[23H25] , nano-mechanical oscillators [53 [59] , Josephson 
junctions [60H62] , and optomechanical crystals [63] . 

Consider a chain of coupled harmonic oscillators with 
Hamiltonian 



N N-l 



i=l 



= 1 



K(a\ai 



+1 



(23) 



As before, we begin by diagonalizing the Hamilto- 
nian. Let us define bk = 3[ sm ^qpr a j> w ^ A — 
y/(N + l)/2 and k = 1, • • • , TV, yielding H = £^(0; + 
£k)b\bk, where = 2ft cos (-^j). The perturbative 
coupling of the two additional quantum registers at the 
ends of the oscillator chain is given by, H' B = g{a\a\ + 

ajy-a/v+i + h.c.) + bj'(a\a$ + a} N , 1 ajv+i), where g charac- 
terizes the register-oscillator-chain coupling strength and 
uj' is the register frequency. Upon re-expressing a\ and 
aiv as a function of the eigenmodes we arrive at the 
full Hamiltonian, 

N 

H B + ^ = ^tfc(46 fc + (-l) fe " 1 aj Nr+1 6 fc +h.c.) 



N 



) + ^(a; + € fc )6l6fc, (24) 



k=i 



where we let tk = (g/A) sin[/c7r/(TV + 1)]. In analogy to 
Sec. [H] we consider resonant tunneling through a partic- 
ular mode b Zl by tuning uo r = uo + e z and ensuring that 
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t z <C \ e z — e z ±i|- The resulting effective Hamiltonian is 
Hf f} = \/2t z (rilb z +blrio), where t? = l/v^(ao + ajv+i)- 
To demonstrate state transfer, we introduce operators 
£± = l/v^^o ± 6 Z ), yielding 



#5/ 



(25) 



Let us now consider unitary evolution under for 

a time r B = ir/(y/2t z ), wherein U^ ff = e~ iH -ff TB = 

so that (U* f )H±(UZ f ) = Re- 
turning to the original basis and evaluating the time evo- 
lution of ao and a/v+i yields 



a (r) - 
«at+i(t) 



(/7f // )ta (/7f // ) = -a 7V+1 , 
^(/7f // )W +1 (/7f // ) = -a , 



(26) 



demonstrating a swap gate between the oscillator- 
registers at the ends of the chain. As before, this state 
transfer is achieved independent of the state of the in- 
termediate chain. Moreover, there exists no additional 
entanglement between the registers and the intermediary 
oscillators; this is a direct consequence of the bosonic na- 
ture of the modes, which, unlike their Wigner-fermionic 
counterparts in Sec. [TTJ carry no strings. 

One crucial difference with the spin-chain case is that 
the occupation of the bosonic eigenmodes is not limited 
to or 1. In a highly excited system, this induces a 
"bosonic enhancement" of off-resonant errors and will 
limit the achievable state transfer fidelity as a function 
of temperature. In particular, the state transfer uni- 
tary evolution gives cln+i(t) = Mjv+i,o a o + V^ a o where 
e = 1 — |Mat+i,o| 2 oc g 2 is a small error and a e is a normal- 
ized linear combination of the modes (i = 1, . . . , N+l). 
The total number of excitations in mode N+l after the 
state transfer is (njv+iM) = (1 — e)(no) + c(n e ), where 
Hi = a\di. Therefore, if the chain is thermally occu- 
pied with (n e ) « kT/uj > 1, the coupling strength g 
must be reduced to g^/uo/(kT) in order to keep errors 
comparable with the zero-temperature bosonic case. In 
realistic experimental systems, this implies an interplay 
between temperature, which sets the bose-enhancement 
of off-resonant errors and decoherence rates, which limit 
the minimal speed of state transfer. 
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FIG. 3: (color online), (a) Contour plots for N = 11 charac- 
terizing the average achievable fidelity as a function of the NV 
depolarization time (Ti) and the coupling strength disorder 
induced by imperfect implantation. Numerics utilize an aver- 
age intrachain spin spacing of d — lOnm corresponding to a 
k — 50kHz dipole-dipole interaction strength. This intrachain 
spacing is assumed to be independent Gaussian distributed 
and the implantation deviation represents the standard de- 
viation, ad- For each cr^, 1000 realizations were averaged to 
obtain the plotted fidelity and a smooth contour plot is gener- 
ated via a third-order spline interpolation. The register-chain 
coupling strengths gi,, gR, Eq. (29), are assumed to be fully 



tunable via control of the 3-level NV ground state manifold 
[27] . (b) Analogous contour plots for N = 51. In this case, 
the NV registers are separated by order optical wavelength 
enabling individual laser manipulation without the need for 
subwavelength techniques, (c) Participation ratio for eigen- 
modes (N = 51). Each eigenmode is indexed by its PR and 
the number of states within a certain PR bin is shown. For 
each disorder (which are represented as fractions of the bare 
coupling strength k = 50kHz), 1000 realizations are averaged. 



Eigenmode-mediated state transfer naturally finds use 
in a variety of quantum computing architectures where 
data buses are required to connect high-fidelity remote 
registers [3j |U l45l . Within such architectures, it is cru- 
cial to consider an interplay between naturally occurring 
disorder and finite decoherence rates. While disorder in 
ID systems generically localizes all eigenmodes, leading 
to an exponentially long state-transfer time, in finite-size 
systems with weak disorder, the localization lengths can 
be large relative to the inter-register separation. In these 



cases, one must still reduce the register-chain coupling 
strength g to compensate disorder effects, but so long as 
the register decay time is sufficiently long, it remains pos- 
sible to achieve high-fidelity QST. In this section, we will 
discuss the impact of coupling-strength disorder on spin 
chains and will analyze the optimization of g as a func- 
tion of disorder strength and qubit depolarization time. 

In particular, we will consider two sources of error: 
1) off-resonant coupling to alternate eigenmodes (which 
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becomes enhanced as disorder increases) and 2) a finite 
register depolarization time Ti, 



2 I^mJ 



2 l^fc.fl 



2 t 



(27) 



where ^(i?) is left (right) register-chain coupling, ipk,L(R) 
is the eigenmode amplitude at the left (right) register, 
Afc is the energy difference from mode z to mode /c, TV 
is the chain length, t is the state transfer time and T\ is 
the depolarization time of the register. The additional 
factor of N in the final term results from the entangle- 
ment discussed in Sec II; indeed, since each register is 
CP-entangled with all intermediate spins, any spin-flip of 
the intermediate chain immediately dephases the quan- 
tum information. 

To ensure that the tunneling rates at each end of the 
intermediate chain are equivalent, we envision tuning 
gL and gn independently, such that t z = gh\^z,h\ — 
9r\^z,r\- Plugging in for the state transfer time, t = 
7r/y/2t z yields, 



k^z 



I^.lI 2 



I^,l| 2 \rl>k,n\' 
\^,R\ 2 Af 



Ntt 



^t i9l \^ l \ j 

(28) 

which enables us to derive the optimal coupling strength, 



9l 




\ 2 x /2T 1 |^,l 



1. Disorder Numerics for a Specific NV-based Architecture 

We now consider an example implementation of 
eigenmode-mediated state transfer in the context of a 
quantum computing architecture based upon Nitrogen- 
Vacancy (NV) registers in diamond [47H49] . Each fully 
controllable NV register consists of a coupled electronic 
and nuclear spin. The nuclear spin, with extremely long 
multi-second room-temperature coherence times is often 
thought of as the memory qubit [8], while the electronic 
spin, which can be optically initialized and read out, me- 
diates interactions with other NVs [3j [45] . Our analysis 
of disorder effects will be based upon the specific archi- 
tecture proposed in [3j; there, NV registers are connected 
by a dark- spin- chain data bus composed of spin- 1/2 elec- 
tronic spins associated with Nitrogen impurities. One of 
the crucial advantages of utilizing spin chains to connect 
remote NVs is that this enables optical addressing of in- 
dividual registers in parallel, a necessary requirement for 
scalable fault-tolerant quantum computation. 

We consider realistic experimental parameters, with 
an average spin spacing of about lOnm, correspond- 
ing to a dipole coupling strength & 50kHz. At room- 
temperature, NV centers are characterized by T\ ~ 10ms 
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FIG. 4: (color online). Coupling pattern {Ji} between spins 
for two differing cases: 1) engineered couplings (circles) as 
in [31] and 2) strong coupling regime (squares). The left y- 
axis characterizes the coupling strength for each case and is 
associated with solid symbols; the couplings are plotted be- 
tween spin numbers (e.g. Jo is plotted between spin number 
and 1). The right y-axis characterizes the fermionic spectrum 
(in this case, the x-axis is simply an index) and is associ- 
ated with the open symbols. The open red circles depict the 
exactly linear spectrum of engineered chain, while the open 
green squares depict the quasi-linear spectrum of the strong 
coupling case with uniform interchain couplings k = 1 and 
optimized g ~ 0.7. 



[8], owing to an Orbach spin-lattice relaxation process; 
the exponential dependence of the Orbach process on 
temperature suggests that slight cooling can significantly 
extend Ti, with many seconds already demonstrated 
at liquid Nitrogen temperatures [3] [68] . We now per- 
form disorder- averaged numerics for two separate chain 
lengths: 1) sub-wavelength addressable (N = 11) and 2) 
optical- wavelength addressable (N = 51) [64]. We char- 
acterize the amount of disorder by the standard deviation 
associated with imperfect spin positioning; in the case of 
NVs, the origin of this imperfection is straggle during the 
ion-implantation process [65, 66 . We average over 1000 
disorder realizations and calculate the fidelity, 1 — 6, ac- 
cording to Eq. (28); in particular, for each realization, we 



calculate the error for each eigenmode of the spin-chain 
and choose the maximum achievable fidelity. As shown 
in Fig. [3^i, high-fidelity quantum gates can be achieved 
for few nanometer straggle provided that the NV depo- 
larization time is ~ 200ms; similarly, for the longer chain 
case (Fig. ^jp) with N = 51, high-fidelity gates are also 
possible, but require significantly longer T\ of a few sec- 
onds. 

Next, we analyze the participation ratio (PR), 



N PR 



1 



(30) 



which provides a characterization of the number of sites 
which participate in a given eigenmode; modes are typi- 
cally said to be extended if Npr ~ O(N) and localized if 
Npr <C N. As the disorder increases, Npp drops sharply 
as depicted in the histograms in Fig. |3fc. Moreover, by 
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<j k « 0.5ft, on average, all eigenmodes exhibit a state 
transfer fidelity < 2/3 even for extremely long T\ ~ 5s. 



VI. STRONG REGISTER COUPLING 

The eigenmode-mediated QST discussed above oper- 
ates in the weak coupling regime, gip <C n/N. Numer- 
ical simulations reveal that by optimally tuning g = 
9m{N) ~ ft, high-fidelity QST can also be achieved 
(see Fig. [5]) . This "strong-coupling" regime enables faster 
state transfer and has been discussed in several recent 
studies [I8j [33 [69] [70] , which focus on the case of an ini- 
tially polarized intermediate chain. Here, we will demon- 
strate that chains with infinite spin-temperature can nev- 
ertheless support QST in the strong-coupling regime. 

To provide intuition for this strong-coupling regime, 
we will begin by considering the engineered spin-chain 
described in [31], where we have N + 2 spin- 1/2 
atoms with nearest-neighbor XX-interactions. The intra- 
chain coupling is non-uniform and is given by, Ji = 
\ yJ(i-\- 1)(N + 1 — z), yielding a Hamiltonian 



N 



#=^J;(<7+^ +1 +/l.C.) + 



N+l , 
i=0 



(31) 



where h is a uniform background magnetic field. Upon 
employing the Jordan Wigner transformation, we once 
again return to a simple tight-binding form, with H = 
J2ij K ij4 c j where = JiS jyi+1 + JjSij+i + hS iyj up 
to a constant. Diagonalizing reveals H = Wkflfk 
with a linear spectrum given by ujk = k + h — . 

As described in Sec. |III[ the system's evolution is gov- 
erned by a(t) = J2j Mij(t)ci(0). Upon setting h = 
one finds that at time t = 2ir, M(2ir) = I and there- 
fore Ci(2ir) = q(0), returning the system to its initial 
state. As the coupling pattern {Ji} harbors mirror sym- 
metry with Ji = Jjv-i, the orthogonal transformation, -0, 
which diagonalizes H can also be chosen mirror symmet- 
ric, i/; ik = (-l^+ViV+i-^. Setting h = |(iV+ 1) 
and t = 7r yields, 



Ma, 



5> 



N+l 



'jk 



(32) 



To demonstrate state transfer, let us recall the analytic 
single- swap fidelity given by Eq. (20). For the moment, 
let us assume that the spins {0, 1, . . . N} are all polarized, 
so that Tr [pf,fP] = 1. Combined with Eq. which 
ensures M ,jv+i = 1, we find Fss = 1, enabling perfect 
QST. We note that in lieu of applying a uniform magnetic 
(iV+1), one can also just apply a simple phase 
1 
(-i) N+ \ 

Turning now to the case of an unpolarized spin chain, 
we again employ the two-qubit encoding previously de- 
scribed. In this case, one will need to apply the phase 



field h=l 
gate Up = 



on spin following transfer. 
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FIG. 5: (color online). Strong coupling regime: By tuning 
g/n ~ iV -1//6 , we obtain high-fidelity QST utilizing an unpo- 
larized chain with two-qubit encoding (paired-protocol). The 
transfer time scales linearly with N (Lieb-Robinson bound) 
[67] and high fidelities > 90% can be maintained for chain 
lengths up to N = 100. 



gate, Up = ^ ^ -^yv+i^ to the logical qubit after state 
transfer. 

The state transfer fidelities for these two strong cou- 
pling methods are given analogously by, 



Fss 



1 1 



[2|M ,iv+i| + |Mo,iv+i| 2 ), 



(33) 



2 6 

Fenc = ^ + ^[2|M , iV +i| 2 |M 2 iV+1 -Mo,oM i v + i,iv+i| 

N 

| 2 + l5^Mjv+i > iM < , | 2 ]. (34) 

i=l 

While these expressions are valid for an arbitrary cou- 
pling pattern (so long as the resultant fermionic Hamilto- 
nian is quadratic), to ensure high-fidelity QST, we require 
|^o,at+i| ~ 1- As depicted in Eq. ( |32| ), satisfying this 
constraint is intimately related to the linear spectrum 
resulting from the choice of Ji = \ \J\i + 1)(N + 1 — i). 

Let us now consider the strong coupling regime (g ~ 
k) where Jo = Jn = 9 and J\ = J 2 = • •• = Jn 1 = 
k. Surprisingly, tuning only g/n enables one to obtain 
a quasi-linear spectrum [18]; such a spectrum will then 
ensure that |Mo,at+i| ~ 1, as desired. Of course, for N = 
2,3, Ji = \y{i + l)(iV + 1 — i) can be satisfied exactly. 
Although for N > 3, an exactly linear spectrum cannot 
be obtained, it is possible to optimally tune g = gM(N), 
so that ujk looks nearly identical to the previous linear 
spectrum, k — (h = 0), as shown in Fig. [I] In 

particular, by optimizing F enc: we obtain gM ~ A/" -1 / 6 , 
with a transfer time r ~ N (Fig. [5|, consistent with [18]. 
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FIG. 6: (color online). Infidelity of QST for a strongly 
coupled chain with long range interactions (1/Ti is assumed 
negligible). Encoded state transfer (paired protocol) fideli- 
ties are shown for dipolar (diamonds), NNN-canceled-dipolar 
(squares) and NN interaction (circles) models. 



VII. LONG-RANGE INTERACTIONS 

Finally, we now consider the effect of longer range in- 
teractions. The majority of proposals for spin-chain state 
transfer focus on approximate nearest-neighbor models; 
however, the microscopic magnetic dipolar interaction is 
naturally long-range and decays as 1 /r 3 , inducing an im- 
portant infidelity in quantum state transfer. The ori- 
gin of this infidelity becomes especially evident as we 
examine the Jordan- Wigner fermionization of the spin 
chain. Each Wigner fermion carries a string of the 
form e l7T a j a j . In the nearest-neighbor case, all 
such strings cancel pairwise leaving a simple quadratic 
model. However, with longer-range interactions, un- 
canceled strings remain and generically introduce pertur- 
bative quartic terms into the Hamiltonian. These quar- 
tic terms imply that the model, unlike the transverse 
field Ising model, is no longer diagonalizable in terms of 
free fermions. In the previous free fermion case, the en- 
ergy of each eigenmode is independent of the occupation 
of all other eigenmodes; this enables state transfer even 
when the spin-temperature of the chain is effectively in- 
finite. By contrast, the quartic terms associated with 
the long-range dipolar coupling introduce interactions 
between fermionic eigenmodes; the energy fluctuations 
of each eigenmode, caused by changing occupations of 
other modes, naturally dephases quantum information, 
limiting the operational spin temperature of the chain. 

Certain proposals have suggested the possibility of us- 
ing dynamical decoupling to effectively cancel next-to- 
neareast neighbor (NNN) interactions [3], but the com- 
plete canceling of all long-range interactions requires a 
level of quantum control that is currently beyond the 
realm of experimental accessibility. Since any long- 



range XX coupling destroys the quadratic nature of 
the fermionic Hamiltonian, an analytic solution for state 
transfer fidelities in the presence of full dipolar interac- 
tions is not available. Thus, we perform exact diagonal- 
ization for chains of length up to N = 12 (total number of 
spins), as shown in Fig. [6] We obtain the encoded state 
transfer fidelities for dipolar, NNN-canceled-dipolar and 
NN interaction models. Remarkably, even with full dipo- 
lar interactions, fidelities ~ 90% can be obtained for a 
total of N = 10 spins; in the case where NNN inter- 
actions are dynamically decoupled, the fidelities can be 
further improved to ~ 98% at similar lengths. 



VIII. QUANTUM MIRROR ARCHITECTURE 

In this section, we present an alternate quantum com- 
puting architecture based upon pulsed quantum mirrors 
[40| |4T] . By contrast to eigenmode- mediated state trans- 
fer, remote quantum logic will be achieved by global ro- 
tations and NN Ising interactions. To remain consistent, 
we choose to discuss the advantages and disadvantages of 
such an architecture within the context of NV registers. 
In particular, analogous to Sec.[VJ we consider NV regis- 
ters connected by spin 1/2 chains of implanted Nitrogen 
impurities. 

Let us begin with a detailed discussion of the mixed 
spin system composed of NV centers and Nitrogen im- 
purities [3]. The full Hamiltonian of a single Nitrogen 
impurity is, 



H 



N 



<y e B-S- lN B • / + A\\ S Z I Z + A ± (S X I X + S y I y ) , 

(35) 

where S is the spin-1/2 electronic spin operator, / is 
the nuclear spin operator, and An = —159.7MHz, A± = 
— 113.8MHz are the hyperfine constants associated with 
the Jahn- Teller axis. 

We envision the application of a magnetic field and 
field gradient, which, within a secular approximation, re- 
duces the Hamiltonian of a nearest neighbor Nitrogen- 
impurity chain to Ising form [3], 



N-l 



H N = k ^2 



N 

£ 



(ujo + Si)Sf, 



(36) 



where n is the relevant component of the dipole ten- 
sor, ujq captures the electronic Zeeman energy, and 5{ 
characterizes the hyperfine term, which is nuclear-spin- 
dependent, for each impurity. Taking into account the 
magnetic dipole coupling between the electronic spin of 
the NV register and the surrounding Nitrogen impurities 
allows us to consider the mixed spin system, 

a-1 N-l 

H eff = ]T KSfSi +1 + JS z NV {Sl + s£) + J2 kS?S? +1 , 
1=1 i=b 

(37) 

where J is the strength of register-impurity interac- 
tion, the Zeeman term in Eq. (36) is assumed to be 
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FIG. 7: (color online), (a) In a qubit chain of length TV, a 
mirror swap operation is defined as the pairwise swap between 
the (l,iV), (2, TV — 1), • • • qubits. This can be achieved via 
global control in the form of single qubit gates (Hadamards) 
and controlled phase gates. Regardless of the initial state, a 
mirror swap occurs after a N + 1 cycles of Q — H • CP [40] . 
where H represents a global Hadamard operation and CP 
denotes a global controlled phase operation, (b) Coupling the 
central NV register to the NV on the left (right) requires the 
ability to perform a directed swap to a neighboring Nitrogen 
spin (grey) . These directed swap operations are made possible 
by using combinations of the mirror sequences Qm (swapS a 
pair of impurities directly surrounding an NV) and Ql (mirror 
swaps an individual impurity chain), (c) Utilizing a fast echo- 
pulse on the NV register (in combination with Qm and Ql) 
allows one to generate selective interactions between the NV 
and any outlined Nitrogen. 



echoed out, and superscripts a, b represent the pair 
of nearest-neighbor impurities next to a given register 
(assuming for simplicity a ID geometry as shown in 
Fig. [7]). The selective individual addressing of the NV 
registers is accomplished via a combination of optical 
beams and microwave driving; this enables an isola- 
tion of the coupling between the NV register and the 
two neighboring impurities. In particular, it is possi- 
ble to perform unitary evolution of the form U e ff = 



e -iH eff T'/2 S ^ ve -iH eff T'/2 S . 

hence, 



NV 



e iK ^ s i s i+i T ' and 



Ui 



local 



= e ~iH eff T e -iK^S*S* +1 T' = e -iJS z NV {S z a +S z h )T 




Q : Nitrogen impurity 
^ : Two-qubit NV Register 



(38) 

by choosing k(T + T') = 2i\m for integer m. We note 



FIG. 8: (color online). Schematic diagram of the 2D compu- 
tational lattice showing limitations imposed by missing Nitro- 
gen implantations (stars) and imperfect Nitrogen-to-NV con- 
versions. Coherent coupling of distant NV registers in a faulty 
2D array can be achieved via global pulsed control of a spin- 
chain quantum data bus. A combination of optical beams 
and a magnetic field gradient allows for individual control of 
NV registers; combined with global single qubit operations on 
Nitrogen impurities in any given row (orthogonal to the field 
gradient) [3 , directed swap operations (e.g. dark green path) 
can be achieved, which allow for quantum information trans- 
fer along arbitrary paths. This field gradient enables a swap 
gate to be performed between two NV registers in adjacent 
rows, which occupy the same column. Moreover, it in fact 
also enables any pair of rows to be swapped, provided that 
the intrarow interactions refocus. 



that this condition implies that the fidelity of Ui oca i is 
extremely sensitive to both coupling-strength disorder as 
well as the general long-range nature of the dipolar in- 
teraction. 



A. Globally Controlled Mirror swap 

Considering only global addressing of the Nitrogen spin- 
chain and unitary evolution as described above, we 
demonstrate a universal set of operations between remote 
NV registers. Coherent register coupling is achieved by 
means of global pulses which mirror the quantum state 
of the impurity chain [40 ; the pulses take the form of 
Hadamard gates and controlled phase gates, which can 
be generated by evolution under an Ising Hamiltonian. 
In an impurity spin-chain of length iV, the global pulses 
swap the state of the first and N th spin, the state of 
the second and (TV — l) st spin etc, as shown in Fig. 
[7^l. The total mirror swap results from TV + 1 cycles 
of Hadamard and controlled phase gates on all impuri- 
ties, Q n+ i = (n^rll^) n+1 . This globally controlled 
impurity mirror will ultimately enable the directed and 
coherent interaction between remote NV registers. 

Let us now consider a specific NV register, separated 
from neighboring registers by impurity spin-chains on 
both sides, as shown in Fig. [7}}. Since the Ising Hamil- 
tonian generates a controlled phase gate, it is possible 
to achieve a mirror swap between any set of qubits con- 
nected by Ising interactions. In particular, the impurity 
Ising interaction allows for a mirror operation within any 
impurity chain, while the Ising interaction corresponding 
to U i oca i allows for a three qubit mirror centered around 
any NV register. This local unitary enables an operation, 
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Qm, which swaps the state of the Nitrogen neighbors of 
the central NV register as shown in Fig. [TJd. 

To couple the central NV register to a specific side 
register, it will be necessary to break the left-right sym- 
metry of the Ising interaction; this is achieved by ex- 
ploiting the length asymmetry between Nitrogen chains 
to the left and right of the NV register. Indeed, it is often 
possible to refocus the mirror operation in one impurity 
chain while causing the edge impurity pair to swap in 
the other chain; we will denote this operation as Ql, as 
shown in Fig. [7}}. Combinations of Qm and Ql success- 
fully manipulate and permute the impurities such that 
the nearest neighbors of the central NV register can be 
any pair of the three impurities (blue, red, green), as de- 
picted in Fig. [7)3. In combination with local rotations of 
the central register, this enables the application of "di- 
rected" unitary evolution, e.g. Udirected — e~ USNvSN b T , 
allowing for the NV register to selectively couple to ei- 
ther side. This enables an interaction between any pair of 
neighboring NV registers effectively mediated by a single 
Nitrogen impurity, 

Hmed ~ J(SnV! + ^NV 2 )^N b ^ (39) 

where NVi and NV2 denote the neighboring registers to 
be coupled and represents the mediating impurity. 
The form of this Ising interaction implies that an ap- 
plication of Qm on this effective three qubit system will 
swap the quantum information of the two electronic spins 
of the remote NV registers. Since each NV center har- 
bors a nuclear-spin qubit in addition to its electronic spin 
[47] . the "there-and-back" protocol described in Sec. [ll| 
enables universal logic between remote registers. 

Having achieved the ability to coherently couple dis- 
tant NV registers within a row, assisted by Nitrogen 
impurities, we now turn to the coupling between adja- 
cent rows in a two-dimensional lattice (Fig. [8]). The sim- 
plest approach involves applying a magnetic field gradi- 
ent along the columns. This would enable a swap gate to 
be performed between two NV registers in adjacent rows, 
which occupy the same column, provided all other inter- 
actions are echoed out. The limited occurrence of ver- 
tically adjacent NVs is a significant source of overhead; 
however, this limitation can be overcome if we achieve 
the ability to swap any pair of nearest neighbor qubits 
in the two-dimensional array, essentially allowing for the 
construction of arbitrary paths (Fig. [8J. Moreover, the 
ability to swap along arbitrary paths also provides an el- 
egant solution to the experimental limitation imposed by 
implantation holes, where a Nitrogen impurity may be 
missing from the ideal 2D lattice. Finally, it also enables 
the use of nominally dark Nitrogen impurities as com- 
putational resources, thereby significantly increasing the 
number of effectively usable qubits. 

While arbitrary individual control of impurities would 
trivially enable such a scheme, realistic constraints limit 
us to individual control of NV registers and only global 
control of the impurity chains. Thus, it is necessary to 
utilize the permutation operation inherent to individual 
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FIG. 9: (color online). Individual control of any given NV 
register (row 1) enables a swap operation between any two 
neighboring qubits along the same row. We illustrate the 
specific example of a four qubit chain. The depicted gate se- 
quence achieves a swap gate (up to individual qubit rotations 
A, B) between the second and third qubit by only applying 
a local gate X on the first qubit (NV register) and global 
operations elsewhere. 

cycles (n^i ' nCi^) n+1 of the mirror operation. These 
gate cycles correspond to an effective propagation of lo- 
cal gates via a relabelling of qubits within a given chain. 
In the simplest scenario, it is possible to apply a swap 
gate between the second and third qubit by only utiliz- 
ing local rotations on the first qubit and global operations 
elsewhere, as shown in Fig. [9] The fundamental opera- 
tion to be propagated is U p = CP -X\-CP where X\ is an 
x rotation (by n) on the first qubit and CP represents a 
global controlled phase gate; propagation takes the form 
of conjugation by mirror cycles where Qk = (H • CP) k 

(k) t 

and Up = Q k U p Qk- To apply a swap operation on the 
n and n + 1 qubit, we let k = n — 1 and apply 

U swap = HU^HXU^ZHU^H, (40) 

where X is a global x rotation and Z is a global z rota- 
tion (by 7r). This protocol requires the ability to produce 
a boundary at the location of the first qubit and allows 
for swaps between arbitrary spins in a given row; moving 
quantum information between rows can be achieved pro- 
vided intrarow interactions refocus (e.g. if vertical and 
horizontal nearest-neighbor distances differ). 

IX. CONCLUSIONS 

In summary, we have extended the analysis of 
eigenmode-mediated state transfer to a variety of im- 
perfections ranging from disorder-driven localization to 
uncompensated long-range interactions. By calculating 
the analytic channel fidelity associated with eigenmode- 
mediated state transfer, we clarify the effects of entangle- 
ment arising from the protocol and illustrate the method 
in which the two-qubit encoding overcomes this chal- 
lenge. We analyze our protocol in the context of pro- 
posed solid-state quantum computing architectures; nu- 
merical simulations with realistic experimental parame- 
ters reveal that QST errors can be kept below certain 
surface-code error-correcting thresholds. Furthermore, 
we have generalized our protocol to the case of bosonic 
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oscillator systems. This approach may enable the rout- 
ing of a "ground-state-cooled" mode through a relatively 
"hot" intermediate oscillator chain, thereby significantly 
reducing the resources associated with system- wide cool- 
ing. 

Moreover, our work may also provide insight into gen- 
eralized infinite-temperature state transfer. In partic- 
ular, by introducing a time-dependent control of the 
register-chain coupling, one may be able to compensate 
for off-resonant errors. This approach finds analogy to 
the continuum wave-packet limit, where dispersion limits 
transfer fidelities; in this case, pre-shaping of the packet 
can overcome nonlinearities of the dispersion. 

Finally, we describe an alternate architecture based 
upon global control pulses which also enables remote 
quantum logic; in particular, we demonstrate that even 
intermediate chain spins can be used as registers, de- 
spite the fact that they are unable to be individually ad- 
dressed. This may provide the blueprint for a novel quan- 
tum computing architecture which utilizes dark spins as 
quantum memory resources. 



X. ACKNOWLEDGEMENTS 

We gratefully acknowledge the insights of and con- 
versations with Pierre Meystre, Alex Zhai, Sidd 
Viswanathan, Peter Komar and Hendrik Weimer. 
This work was supported by the NSF, DOE (FG02- 
97ER25308), CUA, DARPA, AFOSR MURI and NIST. 
AVG acknowledges support from the Lee A. DuBridge 
Foundation. LJ acknowledges support from the Sher- 
man Fairchild Foundation and the NBRPC. AVG and 
LJ acknowledge the IQIM, an NSF Physics Fron- 
tiers Center with support of the Gordon and Betty 
Moore Foundation. The Michigan group was sup- 
ported by the ARO and the AFOSR MURI program, 
the IARPA MUSIQC program, the DARPA OLE pro- 
gram, and the NBRPC(973 Program) 2011CBA00300 
(2011CBA00302). 



Appendix A: Perturbative Calculation of Channel 
Fidelity 

As an extension of the analytic fidelity derivations pre- 



sented in Sec. |III[ here, we provide a closed form ex- 
pression for certain relevant matrix elements. We will 
work perturbatively (g <C n/y/N) and will begin with 
the case of odd chain length. Let us consider comput- 
ing 1 - |(0|M|JV + 1)| 2 for z = (N + l)/2. Recall that 
K is the (N + 2) x (N + 2) coupling matrix of the full 
XX Hamiltonian. We can characterize it with basis \j) 
(j = 0, 1, N + 1) and express 

N-l 

K = 5 (|0)(1| + \N)(N + l\ +h.c.) + + +h.c), 

(Al) 



where we have set the intrachain coupling strength n = 1 . 
The time required for eigenmode-mediated state trans- 
fer is t = V7V + 171-/(20) and M = exp(-iKt). Let us 
now define |±) = (jo) ± |JV + l))/y/2 and further sup- 
pose that TV = Ax — 3 for x G Z>q (we will consider the 
other case below). After going into the diagonal basis \k) 
(k = 1, • • • ,7V) of the intermediate chain, the Hamilto- 
nian breaks down into two decoupled parts as follows: 

K = K + +K_, (A2) 

K_ = [A fc |fc)(fe|+f2 fc (|-)(fc|+h.c.)], (A3) 

fc=even 

k+ = J2 [^k\k)(k\ + n k (\-)(k\ + h.c.)], (A4) 



fc=odd 



where At 



2cos[7r/c/(iV + 1)] and Q k 



(2g/^N + l)sin[7rfe/(iV + 1)]. First consider 
which deals only with even k and does not deal with the 
zero-energy mode \z) = |(iV + l)/2). The eigenstates are 
perturbed only slightly from the original states and we 
call them |— ) and \k) with energy and A& + 0(g 2 ), 
respectively. Moreover, we have 



I-) 



So 



He" 



iK-t 




: >-E^>- ( A5 ) 



2 E(^) (1-cos(A,«)), (A6) 

k<z v k/ 



where as discussed, the sum here is only over even k. 

We now consider K+, which deals with odd k and 
is a little more difficult to treat since it has the zero- 
energy mode. The eigenstates are \s) (for symmetric) 
with energy f] + 0(g 3 ), \a) (for antisymmetric) with en- 
ergy — Q + 0(g 3 ), and |fc) (for all odd k except for k = z) 
with energy A& + 0(g 2 ). We find 

w-£(-E(£)> + »)-g£». 

(A7) 



So 



<+|e-**+*|+> « -1+2 £ (^) 2 (l + cos(A fc 



*)), (A8) 



k<z 

where the sum is only over odd k. 

Putting the results together, we obtain 

1 - |(0|M|iV+l)| 2 

= i-i|( + | e -^ + *| + )-(-| e -^-*|-)| 2 



1 " \ (" 2 + 2 E (^) 2 (! " H^ C0S ( A ^)) 

2 Ef^) 2 ( 1 -(- 1 ) fecos ( A ^))' ( A9 ) 
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FIG. 10: (color online). Depicts a comparison between the 
perturbative (dotted) and exact (solid) calculation for the ma- 
trix element Mo,at+i for N = 205. 1 — |Mo,jv+i| 2 is plotted as 
a function of g, the register-chain coupling strength (k = 1). 
As expected, perturbation theory breaks down as g ~ k/^/N. 
The inset depicts a zoomed in region for 0.005 < g < 0.04, 
where the perturbative expression is in excellent agreement 
with the exact calculation. 



where the sum is over both odd and even k less than 
z = (N + l)/2. Generalizing to all odd N, we find 



|(0|M|7V + 



k<z 



i)l 2 

2 



(1 + (_!)*+* cos(A fc *)), (A10) 



where z = ( N + 1)/2 and where the sum is over both odd 
and even fc (Fig. 10). 

Let us now generalize our expression to include the case 
of even N. Keeping A& = 2 cos[7rfc/(7V + 1)] and Qk = 
(2g/^N + l)sin[7r/c/(7V + 1)], we now have z = N/2, 
t = 7r/O z , and states and N + 1 are shifted by A z + 5 
(note that (5 ~ G(g 2 )). We now define A& = A& — A 2 and 
shift all energies by A z , so that |0) and \N + 1) are now 
at energy 5. Since we are interested in the square of the 
matrix element, this overall energy shift will not affect 
our result. We take K± and |±) as before. As before, we 
first consider N = 4x for x G Z>q (the remaining even 
N will be discussed below). 

We calculate the evolution of |+), which is coupled to 
all odd |fc) and, thus, not coupled to \z). We have 



Ef-i*> 



k \i\ e -iA k t 



(All) 



where A = —St + Y\ t 1 ^ To second order, 



<+| e -"M|+> « e"-£ f^V (l-cos(A fc f)), (A12) 



where the sum is over odd fc. We now consider the evo- 
lution of | — ) (which is coupled to all even including 
We find 



iK_t 




where B = (-5t + £ fc £^) /2. Thus, 

He -i*-*|_)„_ e . 



E 



Ofc 



(l + cos(A fc t)), 



(A14) 

where the sum is over even fc. Putting the calculations 
together, we find 



1, 



1 - |(0|M|iV + l)|^ = l--|(+|e-^|+)-(-|e— |->r 

Generalizing to all even TV, we have, 
1 - |(0|M|iV + l)| 2 



i 



where 



B 



(A - B) 2 
± 

3 + (-l) 2 
4 

3 -(-I) 2 



odd fc^z Afc 



-<ft- 



E 

even fc^z 



of 



(A16) 

(A17) 
(A18) 



Thus, by setting 



(A19) 



i-3(_i)*+* n\ 

h 2 ^ 

we obtain A — B = 0, yielding 
1 - |(0|M|7V + 1)| 2 

w Ef^)^ 1 + (- 1 ) fc+afcos ( A **))' ( A2 °) 

which holds for both even and odd N. One should note 
that tuning S only affects small TV, since for larger TV, 
(A - B) 2 /4 is negligible. 
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We now compute Mq 5 o employing the techniques out- 
lined above; moreover, we note that any requisite matrix 
element entering the average channel fidelity formulae 
can be computed in a similar fashion. For the_Mn o case, 
the K_ expression remains identical to Eq. (A6). How- 



ever, the K+ expression, Eq. (A8), is now different, 

(+\e~ iK+t \+) « l-2g (^j (l-cos(A fe t)), (A21) 

where the sum is over odd k. Combining yields, 
1 



M ,o 



1 - | [<+|e- 



iK+t 



\+) + (-\e 



iK-t\ 



k<z 



(l-cos(A fc t)), 



(A22) 



where the sum is now over both odd and even k. 



Appendix B: Channel Fidelity for Remote <j z Gate 

In this appendix, we illustrate the channel fidelity as- 
sociated with an eigenmode-mediated remote a z gate. In 
combination with the detailed discussion of the double- 
swap channel fidelity in Sec. |III[ this provides the 
framework for calculating the gate fidelity of a remote 
controlled-phase gate. In particular, we examine the 
process whereby: 1) register-0 is swapped across the in- 
termediate chain, 2) a <j z gate is performed at register- 
iV + l and 3) a second return step of eigenmode-mediated 
state transfer is performed. In the ideal case, this remote 
<j z channel should result in a a^-gate on register-0 and 
hence, the associated fidelity is given by 



where we have used the fact that the number of excita- 
tions in each mode must be preserved. As before, for 
i = x,y, only cross terms involving <j + and a~ provide a 
non-zero contribution to Eq. (Bl). We find 



Tr [a z a+a z Ti A [U z (<Tq ® Pch )U}]] 



(M > +1 ) 2 - (M* ) 2 



(MS* 



(B3) 



An analogous calculation yields the contribution from the 
opposite cross term and thus, we now turn to the a z 
contribution. Again, we begin by calculating the time 
evolution, 

<rg(t) = U t e i7vnN + 1 U\2clc -l)Ue i7vnN + 1 U 

-> 2U t e innN + 1 ^M^ i M 0J 4c j e i7rnN + 1 U(B4) 



where we have dropped the (—1) contribution from the 
first line since it will ultimately trace to zero. Conjuga- 
tion by e i7TUN + 1 affects J2i,j M o,i M oj4 c j onl y if i or 3 
equals N + 1; in these cases, the matrix element gets an 
additional negative sign. We can capture this by defining 
an (TV + 2) x (N + 2) diagonal matrix 5, which contains 
unity along all diagonal entries except the last, where it 
contains (—1). Using S, we find, 



crg(t) = 2U^Y J MS :i Moj4c j U 



— » 



2 J2 M^Moj J2 M^M^cUj,, (B5) 

hj i'J' 



F x = 



1 

12 



Y Tr [a z a l a z TT A [U z (a l <g> p ch )Ul]\Bl) where M = MS. A non-zero contribution arises only if 
_ ' ' i' = j' = 0, wherein we find Trfo-Q <7 erg] = 1. Combining 



z=x,y,z 



all contributions yields, 



where U z = U(j z NJrl U, U represents an eigenmode- 
mediated swap, p c h is the mixed state of spins 1, • • • , N+ 
1, Tro traces over register-0 and Tr^ traces over all other j? _ 1 _|_ 1 
spins. Let us begin by calculating the time evolution of 



°o(t) = u t c l u z = U t e i7rnN + 1 U t clUe i7rnN + 1 U 



17+ 



M * o4 - MX 



o,at+i u at+i 



u 



2 6 
1 
6 

ho 

1 1 

2 + 6 



(M > +1 ) 2 - (M * ) 2 - E(M *J 2 + c.c. 
[|(0|MSM|0)| 2 - 2Re((0|MSM |0»] , (B6) 



— » 



(Mq ) 2 Cq — (MQ ?Ar+1 ) 2 cj v+1 + y^(MQ^) 2 cj,(B2) where we have made use of the fact that M is symmetric. 
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